Characteristics of First-Order Vortex Lattice Melting: 
Jumps in Entropy and Magnetization 
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We derive expressions for the jumps in entropy and magnetization characterizing the first-order 
melting transition of a flux line lattice. In our analysis we account for the temperature dependence 
of the Landau parameters and make use of the proper shape of the melting line as determined by 
the relative importance of electromagnetic and Josephson interactions. The results agree well with 
experiments on anisotropic YiBa2Cu307_^ and layered Bi2Sr2CaiCu20s materials and reaffirm 
the validity of the London model. 



PACS numbers: 74.60.Ec, 74.60. Ge 

A cornerstone of the phenomenology of type II su- 
perconductors is the Abrikosov mean-field H - T phase 
diagram comprising a Meissner-Ochsenfcld phase at low 
fields, a Shubnikov- or mixed phase at intermediate fields, 
and a normal metallic phase at high fields H > H C2 
M. Thermal fluctuations may modify this mean- field 
picture considerably as the vortex lattice melts into a 
flux-liquid phase, which is the case in high temperature 
superconductors . According to standard symmetry 
considerations Q, this melting transition is expected to 
be of first order. Recently, this expectation has been con- 
firmed through the experimental observation of a jump 
in the magnetization in layered Bi2Sr2CaiCu2 0s (BiS- 
CCO) |,| and in anisotropic YiBa 2 Cu 3 7 _«5 (YBCO) 
crystals. The latent heat released in the transition 
has been determined in calorimetric measurements on an 
YBCO crystal |||| and the thermodynamic consistency 
between the magnetic and calorimetric experiments via 
the Clapeyron relation has been demonstrated [||. 

Though consistency could be achieved in the experi- 
ments, simple estimates [pj for the jumps in magnetiza- 
tion and entropy have failed to explain the magnitude 
and the temperature dependence of these characteristic 
quantities. In particular, in the layered BiSCCO mate- 
rial the entropy jump per vortex per layer as extracted 
from the magnetization data via the Clapeyron equation 
seems to diverge upon approaching the superconducting 
transition temperature ||, and no explanation of this 
striking result has been given so far. In this letter we re- 
solve this puzzle and derive the temperature dependence 
of the jumps in magnetization and entropy. Our results 
are consistent with all the data measured so far in YBCO 
and BiSCCO single crystals, see Fig. 1. 

Besides construction of simple estimates, the theoreti- 
cal discussion of the jumps in magnetization and entropy 
has concentrated on numerical simulations. The prob- 
lem of the notoriously small entropy jumps obtained in 
early simulations |il]| , |Tl"|] has been cured by the insight of 
Hu and MacDonald p2], who pointed out the relevance 



of the temperature dependence of the Ginzburg-Landau 
parameters in the analysis of numerical data obtained in 
simulations based on the lowest Landau level approxi- 
mation. Their insight can be given a much wider per- 
spective, both in terms of its application to other nu- 
merical approaches Jl3{ as well as analytic treatment of 
the problem: combining this concept with our knowledge 
of the shape of the melting line, we show below how to 
construct a consistent scheme which provides us with ex- 
pressions for the jumps in entropy and magnetization, 
exhibiting both the correct magnitude and temperature 
dependence. Our analysis then removes recent doubts 
|l2| , ^4| on the ability of a simple London model to de- 
scribe the large entropy jumps measured on YBCO and 
BiSCCO crystals close to T c . 

Two important elements in our analysis are the shape 
of the melting line H m (T) and the relevant fluctuation 
mode in the vortex system at the melting transition, the 
latter defining the volume V^f of the elementary degree 
of freedom. Both quantities depend sensitively on the 
degree of anisotropy/layeredness of the material through 
the relative importance of the two types of interactions 
appearing in the vortex system, the Josephson- and the 
electromagnetic interaction [{nil . As a consequence, the 
melting process of the vortex lattice in YBCO and BiS- 
CCO exhibits quite distinct characteristics, both regard- 
ing the shape of the melting line and the temperature 
dependence of the jumps in magnetization and entropy. 
In addition, we will make use of the Clausius-Clapeyron 
equation relating the jumps in the entropy density As 
and in the magnetic induction AB, 



As 



1 dH n 



-AB. 



(1) 



In the following, we first present the analysis for an 
anisotropic material such as YBCO, where we make use 
of the scaling form of the London free energy func- 
tional which contains the flux lattice constant as the only 
length scale, rendering the calculation essentially exact. 
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For strongly layered materials (e.g., BiSCCO), two new 
length scales, the layer separation d and the London pen- 
etration depth A, become relevant. We show how to con- 
struct estimates for the jumps which reduce to the pre- 
vious results in the scaling regime and which allow us to 
deal with layered materials. 

Consider the statistical mechanics of a system de- 
scribed by an effective free energy functional .F(T, </>). 
The coarse grained effective field <j> ( e -g- a Ginzburg- 
Landau wavefunction) is obtained after integration over 
the microscopic degrees of freedom, which generates the 
temperature dependence in T (e.g. the Landau param- 
eter a(T) — a(l — T/T c )). After integration over the 
remaining degrees of freedom <f> we arrive at the parti- 
tion function (k B is the Boltzmann constant) 



Z = V4>e 



-F{T,<t>)/k B T 



(2) 



Taking the derivative of the free energy F = — k B T\nZ 
with respect to T we obtain the entropy, 



S = So - (&rF), 



(3) 



where S = — F)/T is the configurational entropy 
of the coarse grained 0-field, whereas the second term 
accounts for the internal temperature dependence in the 
free energy functional !F(T, (f>). To be specific, let us con- 
sider a vortex system within the London approximation 
where the free energy functional takes the form 

^[{Sf,}} = E£ a / ds M ■ ds v : : — , (4) 



with e = ($ /47tA) 2 the basic energy scale in the prob- 
lem (proportional to the vortex line energy), A the Lon- 
don penetration depth, e 2 — m/M < 1 is the anisotropy 
parameter, and $ = hc/2e denotes the flux quantum. 
In (|) all lengths are measured in units of the vortex 
separation a Q = y/ Q /B, with the configurations ex- 
pressed through the dimensionless position variables s M . 
For the vortex system, the parameters e and A depend 
on the temperature via A 2 (T) = Aq/[1 - (T/T c ) 2 ]. In the 
limit A > a Q the functional (^) assumes a simple scal- 
ing form, with all physical parameters appearing in the 
prefactor se (T)a , the remaining factor representing a 
scale independent summation over geometrical configura- 
tions of lines Jl6|. Within the scaling regime, all physical 
results depend on the combination ee (T)a . E.g., the 
shape B m (T) of the melting line is given by the condi- 
tion e£ {T m )a /k B T m = const. (= 1/2^/ttcI, we make 
use of the usual definition of the Lindemann number c L ), 
from which one easily derives the standard result |5|,|l5|| 



l[ ' A 2 L (fc B T) 2 



Tl) 



(5) 



(note that Eq. (|J) derives from scaling and proves the 
validity of the Lindemann criterion). Returning to (|^), 
we can again make use of the scaling form of T and ex- 
press the second term in (||) through the energy (T): 
(drF) — {dr lne )(jF) and inserting back into (||) we 
arrive at the relation 



5' 



e dT 



T ds F 
~ ~dTT' 



(6) 



At the phase transition the entropy exhibits a jump while 
the free energy F remains continuous. From (^|) we thus 
infer the following relation between the entropy jump AS 
and its configurational part A5 Q , 



AS 
ASo 



T ds 
~ T ~dT 



1 + (T m /T c ) 2 
1 - (T m /T c ) 2 



(7) 



where we have made use of the specific temperature de- 
pendence of the line energy e oc [1 — (T/T c ) 2 ] in the last 
equation. We find, that close to the thermodynamic tran- 
sition the entropy jump AS is strongly enhanced with 
respect to its configurational component AS'o. Note that 
it is the latter quantity which is amenable to simple es- 
timates [[| and which is usually calculated in numerical 
simulations |l3| . The result (Q) can be easily understood 
as the simple consequence of a "coordinate" transforma- 
tion: the temperature enters the partition function Z not 
merely as a scale parameter but rather in a combination 
T/[l — (T/T c ) 2 ]. Close to T c this expression becomes sin- 
gular, resulting in a strong enhancement of the entropy 
jump |TjJ]. The physical origins of this enhanced entropy 
are microscopic fluctuations: Within the coarse grained 
vortex model these fluctuations surface in the tempera- 
ture dependence of the phenomenological parameters in 
T, which are singular at the mean-field transition. 

Next, let us find an expression for the configurational 
part AS'o of the entropy jump: within the scaling regime, 
standard arguments dictate the form 



V 

AS = f)k B ——, 

Vcdf 



(8) 



with Vodf = sal, V denotes the system volume, and r\ is a 
small number. Physically, this result can be understood 
by attributing the thermal energy k B T to each individual 
degree of freedom. The volume 14df is determined by the 
dominant modes leading to melting, which are located at 
the Brillouin zone boundary with k± ps \^7r/a . Hence, 
we can define the volume per degree of freedom in the 
form V c di = a 2 L. Furthermore, in an anisotropic ma- 
terial the important fluctuations involve the wavevector 
k z ~ l/ea along the field, thus L — ea (with all numer- 
icals absorbed in 77). Combining (Q) and (||) we obtain 
the final result for the entropy jump per vortex per layer, 



AS d w 2r? 



ea 1 - {T m /T c ) 2 



(9) 
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On the melting line, the product a [l — (T m /T c ) 2 ] is 
(roughly) temperature independent and we find a con- 
stant but material dependent entropy jump per vortex 
per layer, in agreement with experiments on YBCO [||. 

In order to find a numerical result for the entropy 
jump we compare the latent heat per vortex line Li = 
T m AS a 2 with the result of numerical simulations car- 
ried out within the present London formalism |l3| : Mak- 
ing use of (||) and the numerical result L; = 0.015 s we 
find the expression ASd ~ 0.03 e o (0)d/T c and using pa- 
rameters for YBCO (A w 1400 A and d = 12 A) we 
obtain the value ASd ~ 0.4 k B , in good agreement with 
experiment Using the result for the melting line in 
Ref. 13 1, ee (T m )a / k B T m w 11, we obtain a value for 
the parameter 77, 77 w 0.16. 

Let us turn next to the jump AB in the induction. In 
order to make use of the Clapeyron equation (Q) we have 
to determine the slope drB^ of the melting line |fl8| . 
Ignoring for the moment the temperature dependence in 
the free energy (||) (and thus in the parameters of M)) we 
have dxB m = —2B m /T m . The full result which accounts 
for the temperature dependence in e Q reads 



dB m 
dT 



2B 



T m V So dT J 



(10) 



The same factor relating the two slopes appears in equa- 
tion (Q) which expresses the entropy jump AS through 
AS - Thus in the final expression for the jump AB in the 
induction this correction factor drops out and we arrive 
at the simple result (we use V e df = a 2 L) 



AB 



(11) 



where y, = 2ttt] w 1.0. Using L — ea a as appropriate in 
an anisotropic material we arrive at the final result for 
the jump in B, 



AB 



6.10" 



$0 



A 2 (T m )' 



(12) 



Note that in an incompressible (uncharged, e — > and 
A — > 00) system, we correctly find AB — > 0. Rewriting 
U) in the form AB[G] sa (1.5 ■ 10- 6 /e)T ro [K](B ro [G]) 1 /2 
and choosing e = 1/8 we arrive at a good agreement 
with the magnetization data of Schilling et al. || on an 
YBCO single crystal, see Fig. 1 (we have made use of the 
experimentally measured melting line B m (T)). 

So far our analysis has been essentially exact: we have 
exploited the scaling behavior of the London functional in 
the regime a < A and have deduced the single unknown 
parameter r\ from a comparison to a numerical simulation 
of the London model. While this approach is success- 
fully applied to a continuous anisotropic superconductor 
such as YBCO, we have to reconsider the situation for 
strongly layered materials (e.g., BiSCCO) as new length 
scales (d, A) enter the problem. We then can make use of 



an important insight provided by the above derivation, 
namely that the calculation of the jump in B does not 
suffer from the complications associated with the deter- 
mination of the jump in entropy. This is because the en- 
tropy involves a derivative of the free energy with respect 
to temperature, whereas the induction is given by the 
derivative with respect to the magnetic field. The latter 
usually does not show up in the Landau parameters |l9| . 
Indeed, we can arrive at the result ([n]) starting from the 
general thermodynamic relation B = — (47t/U)<9h-G|t, 
where G is the Legendre transform of the free energy 
F, G(T, H) = F(T, B) - BHV/4ir. With the estimate 
G ~ k B TV/V c df and making use of the power law depen- 
dence of U e df on H we obtain 



AB 



[i k B T n 
HV C 



cdf 



(13) 



Inserting the ansatz V c df — a 2 L — <& L/H we imme- 
diately recover the result (|ll]) (with fi replaced by an 
unknown numerical //). Once the jump in the induction 
AB is known, we can make use of the Clausius-Clapeyron 
equation ([!]) and arrive at the result for the jump in the 
entropy. The problematic factor arising from the tem- 
perature derivative of the free energy functional !F, see 
(||) and (|7[), is taken care of by the derivative drB m of 
the melting line, see ([To|). 

We proceed with the analysis of the jumps in magneti- 
zation and entropy for strongly layered superconductors, 
following the above line of thought. To do so we need 
to know the length L of the relevant modes at melting, 
see (0), as well as the shape of the melting curve in a 
layered superconductor, to be used in (|l|). 

In layered BiSCCO the melting line is pushed down 
to low fields B m (T) < B\(T) = $ /A 2 over a large 
portion of the phase diagram. The dominant interac- 
tion in the vortex system is then given by the electro- 
magnetic one. The loosely bound pancake vortices un- 
dergo large thermal fluctuations and dominate the melt- 
ing process, hence L = d. The shape of the melting 
line follows most easily from a Lindemann analysis with 
(u 2 ) = c 2 al ~ Td/ei(k z ~ while making use of the 

dispersive electromagnetic line tension si ~ e /X 2 k 2 [ ^5| , 



$0 eg e d 

a^Ya^t 



T 2 



(14) 



Close to T c the Josephson interaction becomes relevant 
as soon as eA(T) > d: For T > T om ps T c [l - (eA /d) 2 ] 1/2 
the dominant fluctuations at melting are cut off on the 
larger scale L ~ eA. The Lindemann criterion takes the 
form (u 2 ) = c 2 a 2 ~ TeX/ei(k z ~ 1/eA) and we obtain a 
melting line following a (1 - T 2 /T 2 ) 3 / 2 behavior |]|, 



B 



cm, J 



(T) 



$0 7TC 2 ££ A 

A 2 " 4 k B T 



2\ 3/2 



T 2 ) 



(15) 
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Inserting the above results for L into (|Tl]) we obtain the 
jump AB in the induction for a layered material, 



AB 



T Cm < Tyr 



(16) 



Note that it is the temperature dependence of L, which 
goes from ea (T' m ) in a continuous anisotropic supercon- 
ductor to L = d and L = e\(T m ) in a layered mate- 
rial, that leads to the different dependencies in the jump 
AB(T). 

The result ( |l6| ) explains the experimental observations 
of Zeldov et al. Q, see Fig. 1: At low temperatures 
T m < T c — 7 K the jump AB increases linearly with 
temperature. About 7 K before reaching the transi- 
tion, AB drops sharply and vanishes at T c . This be- 
havior is explained in terms of the crossover at T cm , 
where the Josephson coupling between the layers be- 
comes relevant and cuts off the further growth of fluc- 
tuations. With the same value fi' = fi (rs 1) as above, we 
rewrite (|l6|) in the form AB[G] « (0.07/d[A])T m [K] and 

AB[G] « (0.07/£Ao[A])T m [K][l-(!T m /!r c ) 2 ] 1/2 and using 
the parameters d = 15 A, A ~ 2000 A, and e rs 1/400 we 
obtain excellent agreement with the experimental result 
of Zeldov et al. ||. Note that the experimental find- 
ing AB(T — ► T c ) -> (while B remains > AB), com- 
bined with jll|), hints at a divergence L(T — > T c ) — > oo 
and thus a 3D transition as T c is approached. Finally, 
the expressions for the jumps in entropy follow from the 
Clapeyron equation (using the line shapes (|l4|) and (|loj) 



AS d ^{ 



/' 



7T 1 - (T m /T c ) 2 ' 

3^ d k B 

1 4^^7i = ^7^f' 



T ra < T c 



(17) 



Unlike in the anisotropic case, for the layered material 
the entropy jump per vortex per layer diverges on ap- 
proaching the transition at T c , again in agreement with 
the experimental observation p] (note that at low tem- 
peratures the entropy jump in Ref. || vanishes as the 
melting line flattens, possibly due to disorder). 

In summary, our new estimates , (|l2|) and ( |l6| ) , ( |l7| ) 
provide a consistent explanation of the observed charac- 
teristic jumps at the first-order melting transition of the 
vortex crystal in type II superconductors. 
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FIG. 1. Top: the jump AB in the induction versus tem- 
perature T as measured in an YBCO single crystal [8] and 
calculated from the expression (12). The deviations close 
to the transition are possibly due to sample inhomogeneity 
[7]. Bottom: the same for a BiSCCO single crystal [6] us- 
ing the result (16). The drop in AB on approaching T c is 
explained in terms of a temperature dependent cutoff in the 
electromagnetic fluctuations through the Josephson coupling 
at temperatures T > T om . The entropy jump ASd per vor- 
tex per layer as obtained from AB through the Clapeyron 
equation diverges at the transition [6], in agreement with the 
result (17). 
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